arXiv:1502.03090vl [astro-ph.HE] 9 Feb 2015 


33rd International Cosmic Ray Conference, Rio de Janeiro 2013 
The Astroparticle Physics Conference 



The Statistical and Numerical Study of the Longitudinally Asymmetric Distri¬ 
bution of Solar Proton Events Affecting the Earth Environment of 1996-2011 

Hongqing Weixing Wan' 

' Beijing National Observatory of Space Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, Beijing 
100029, China 

^ CAS Key Laboratory of Geospace Environment, Department of Geophysics and Planetary Sciences, University of Science and 
Technology of China, Hefei, Anhui 230026, China 

hqhe @ mail, iggcas. ac. cn 

Abstract: Large solar proton events (SPEs) affect the solar-terrestrial space environment and become a very 
important aspect in space weather research. In this work, we statistically investigate 78 solar proton events of 
1996-2011 and hnd that there exists a longitudinally asymmetric distribution of flare sources of the solar proton 
events observed near 1 AU, namely, with the same longitude separation between magnetic field line footpoint of 
observer and flare sources, the number of the solar proton events originating from sources located at eastern side 
of the nominal magnetic footpoint of observer is much larger than that of the solar proton events originating from 
sources located at western side. A complete model calculation of solar energetic particle (SEP) propagation in 
the three-dimensional Parker interplanetary magnetic field is presented to give a numerical explanation for this 
longitudinally asymmetric distribution phenomenon. We find that the longitudinally asymmetric distribution of 
solar proton events results from the east-west azimuthal asymmetry in the topology of the Parker interplanetary 
magnetic held as well as the effects of perpendicular diffusion on the transport of SEPs in the heliosphere. Our 
results would be valuable in understanding the solar-terrestrial relations and useful in space weather forecasting. 
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1 Introduction 

Solar energetic particles (SEPs), which are charged ener¬ 
getic particles occasionally emitted by the Sun, risk the 
health of astronauts working in space and damage elec¬ 
tronic components on satellites, so they have become a 
very important aspect in affecting solar-terrestrial space 
environment and space weather. Theoretically, SEPs ob¬ 
served in the interplanetary magnetic Held provide fun¬ 
damental information regarding acceleration mechanisms 
and transport processes of charged energetic particles. 
Therefore, the subject has become a focus of astrophysics, 
space physics, and plasma physics. 

Through several decades of investigations in the com¬ 
munity with spacecraft observations and theoretical mod¬ 
eling, signiHcant progresses have been achieved toward a 
better understanding of the transport processes of SEPs. As 
a pioneering work, m provided a focused transport equa¬ 
tion to investigate the modulation of galactic cosmic rays 
and the transport of SEPs. It is very difficult to analyti¬ 
cally solve the multidimensional transport equation; con¬ 
sequently, numerical calculations are usually adopted (e.g., 

EE]). 

Recently, a investigated the effects of particle source 
characteristics on SEP observations at 1 AU and found that 
the perpendicular diffusion has a very important influence 
on the propagation of SEPs in the heliosphere, particularly 
when a spacecraft is not directly connected to the accel¬ 
eration regions either on the Sun or near the CME-driven 
shocks by the interplanetary magnetic Held lines; in such 
cases, the earliest arriving particles can be seen propagat¬ 
ing toward the Sun, having scattered backward at large dis¬ 
tances. 013 presented a direct method to quickly and ex¬ 


plicitly determine the spatially dependent parallel, radial, 
and perpendicular mean free paths of SEPs with adiabatic 
focusing. They reported that when the turbulence strength 
5B/B is sufficiently large, e.g., dB/B > 2.34, the ratio 
A^/A|| would approach or exceed unity. 

In this work, we statistically investigate 78 solar proton 
events (SPEs) of 1996-2011 observed near 1 AU and And 
that with the same longitude separation between magnetic 
Held line footpoint of observer and flare sources, the num¬ 
ber of the solar proton events originating from sources lo¬ 
cated at eastern side of the nominal magnetic footpoint of 
observer is much larger than that of the solar proton events 
originating from sources located at western side. Our sys¬ 
tematical simulations also confirm the longitudinally asym¬ 
metric distribution of the solar proton events observed in 
the interplanetary space including Earth’s orbit. We find 
that the longitudinally asymmetric distribution of solar pro¬ 
ton events results from the east-west azimuthal asymme¬ 
try in the geometry of the Parker interplanetary magnetic 
Held as well as the effects of perpendicular diffusion on the 
transport processes of SEPs in the heliosphere. Our results 
would be valuable in understanding the solar-terrestrial re¬ 
lations and useful in space weather forecasting. 

2 Data Set 

We collect 78 solar proton events affecting the Earth en¬ 
vironment observed near 1 AU in the time period 1996- 
2011 from Space Weather Prediction Center of NOAA 
at http;//www.swpc.noaa.gov/ftpdir/indices/SPE.txt. The 
data set covers an uninterrupted time period of 16 years 
including the whole solar cycle 23 and the rising phase of 
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solar cycle 24. The maximum proton flux of each particle 
event measured by GOES spacecraft at Geosynchronous 
orbit for energies >10 MeV is greater than or equal to 10 
particle flux units (pfu, 1 pfu=l particles/(cm^ — ir — i)). 
Note that the 78 solar proton events are selected from 
the original listing by excluding the events without defi¬ 
nite location on the solar surface, since the explicit data 
of the active regions are indispensable to the investiga¬ 
tions. The solar wind speed used in this work is averaged 
from SOHO or ACE data measured within the time range 
\tmax — 2hours, tmax + 2hours\, where t^ax is the maximum 
time of each SEP event. 


Solor- Proton Events; 1996-2011 



3 Statistical Results 

The heliospheric magnetic field is “frozen-into” the solar 
wind and travels along with it to form the Archimedes 
spirals. The well-known Parker spiral can be expressed as 

0 , = ( 1 ) 

where Q. is the angular rotation rate of the Sun, r is the he¬ 
liocentric distance, and V™ is the solar wind speed. Eor 
each solar proton event, by inputting the solar wind speed 
y*”' into Equation ([T]), we can derive the Parker spiral an¬ 
gle of the nominal magnetic field line connecting the ob¬ 
server at 1 AU with the Sun. That is to say, we can deter¬ 
mine the heliographic longitude of the nominal magnetic 
footpoint of the spacecraft located at Earth’s orbit during 
the SEP events. Eor convenience in the investigations, we 
mark the west and east heliographic longitudes (j) of the 
flare sources on the solar disk as positive and negative val¬ 
ues, respectively. Then for each solar proton event, we can 
obtain the relatively longitudinal distance 0;- between the 
heliographic longitudes of the flare location and the nom¬ 
inal magnetic footpoint of the spacecraft located at 1 AU 
during the particle event, through the relation 

0, = (^)-(^), = 0-£2r/y™. (2) 

Eigure[T]presents the spatial distribution on the solar sur¬ 
face of the flare sources of the 78 solar proton events ob¬ 
served near 1 AU in the time period 1996-2011. Similar 
to the heliographic longitude, the heliographic latitude is 
marked as a positive value when the flare source located at 
the northern side of the solar equator and a negative value 
when the flare source located at the southern side. Erom 
the upper panel of Eigure [T] we can see that the majority 
of flare sources of the solar proton events are distributed 
in the heliographic latitude range [530°,A30°]. The lower 
panel of Eigure [T] shows the number percentage distribu¬ 
tion of the solar events along the relatively longitudinal dis¬ 
tance which is divided into 18 bins. The event number 
percentage P, (/ = 1,2,..., 18) within each bin of 0,- is ob¬ 
tained as follows: 


where A,- (; = 1,2,..., 18) denotes the corresponding so¬ 
lar event number in each bin. The vertical dashed line in 
the middle of the lower panel in Eigure [T] indicates the 0° 
relatively longitudinal distance, where the nominal mag¬ 
netic footpoint of the observer and the flare source have 
the same heliographic longitude. At the left side of this 
line, the negative values indicate that the flare sources lo¬ 
cate at the eastern side of the nominal magnetic footpoint 


Figure 1: The spatial distribution on the solar surface of 
the flare sources of the 78 solar proton events observed 
near 1 AU in the time period 1996-2011. Upper panel: the 
latitudinal and longitudinal distribution of the solar sources 
of the proton events. Lower panel: the number percentage 
distribution of the solar events along the relatively longi¬ 
tudinal distance (j)r, which is divided into 18 bins. With 
the same longitude separation between the magnetic foot¬ 
point of observer and solar sources, the number of the solar 
proton events originating from sources located at eastern 
(left) side of the magnetic footpoint of observer is much 
larger than that of the solar proton events originating from 
sources located at western (right) side. 


of observer, whereas the positive values at the right side 
of the middle line indicate that the flare sources locate at 
the western side of the magnetic footpoint of observer. It 
can be seen that the event number percentage is maximum 
in the interval [—20°,0°] among all the 18 relatively lon¬ 
gitudinal bins. Interestingly, we find that the event num¬ 
ber percentage distribution is asymmetric about the verti¬ 
cal dashed line in the middle. We connect the event number 
percentages at the central value within each relatively lon¬ 
gitudinal interval with red and blue lines on the left-hand 
and right-hand halves, respectively. We further mirror the 
event number percentages on the right-hand half to the left- 
hand half with a blue dotted line. Then we can clearly see 
that the event number percentage within each interval on 
the left-hand half is systematically larger than that in the 
symmetrical interval on the right-hand half. That is to say, 
with the same longitude separation between magnetic field 
line footpoint of observer and flare sources, the number of 
the solar proton events originating from sources located at 
eastern side of the nominal magnetic footpoint of observer 
is much larger than that of the solar proton events originat¬ 
ing from sources located at western side. 


4 Numerical Model and Results 
4.1 Model 

In our model, the five-dimensional focussed transport equa¬ 
tion that governs the gyrophase-averaged distribution func¬ 
tion /(x,/j,p,f) of SEPs can be written as (e.g., IJl) 
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where x is particle’s position, z is the coordinate along the 
magnetic field spiral, p is particle’s momentum, p is pitch- 
angle cosine, and Q is source term. 

Under the diffusion approximation for a nearly isotropic 
pitch-angle distribution, the parallel mean free path A|| 
could be written as (Il6l) 


All = 


f/ 


3v /■+1 (1 -ju2)2 


D 


-dll. 




(5) 


In this work, we typically use radial mean free path = 
0.3 AU (corresponding to parallel mean free path A|| = 
0.67 AU at 1 AU) and perpendicular mean free paths 
Xx = Xy = 0.006 AU for 50 MeV protons. We utilize the 
time-backward Markov stochastic process method (||2l) to 
numerically solve the focused transport equation (|4|i. In 
our simulations, the Parker interplanetary magnetic field B 
is set so that its magnitude at lAU is 5nT, and the solar 
wind speed is typically set as = 400 km Typical 
sizes of SEP sources (flares or CMEs) are tens of degrees 
wide in heliographic latitude and longitude, so in the sim¬ 
ulations, we set SEP sources with limited coverage of 30° 
in latitude and longitude. Generally, the unit of omnidirec¬ 
tional flux is used as in this 

work, however, we use an arbitrary unit for convenience in 
plotting figures. 


4.2 Numerical Results 

We simulate 18 cases to show the effect of solar source lo¬ 
cation on the SEP flux observed in the interplanetary space. 
All the solar particle sources have the same coverage, but 
their centers are located at various heliographic longitudes 
in the solar equator. Specifically, the longitude separations 
between the centers of 18 solar sources and the magnetic 
footpoint of the observer are set to be the following values: 
0°, 20° west, 20° east, ..., 160° west, 160° east, 180° east. 
All the other conditions and parameters are the same for 
each case. In these simulations, we typically investigate 50 
MeV protons detected at 1 AU, 90° colatitude. 

Eigurel^shows the flux maxima of the SEP events origi¬ 
nating from the 18 solar sources with different longitudinal 
locations relative to the magnetic footpoint of the observer 
at 1 AU, 90° colatitude. The X-axis in Figure |2] is the rel¬ 
atively longitudinal distance (j)r between the solar source 
center and the magnetic footpoint of the observer, same as 
that in Figure [T] The vertical dashed line in the middle of 
Figure|2]denotes the 0° relatively longitudinal distance. As 
we can see, the maximum flux of SEPs in this case is the 
largest among all the 18 cases investigated. This means that 
when the observer in the interplanetary space is connected 
directly to the solar source by heliospheric magnetic field 
lines, the SEP flux is larger than that otherwise. We con¬ 
nect the flux maxima of SEPs in the simulated cases on the 
left-hand and right-hand halves with red and blue curves, 
respectively. It can be obviously seen from Figure |2] that, 
the farther the solar source center is away from the mag¬ 
netic footpoint of the observer, the smaller is the maximum 
flux of SEPs observed (see also Q). Similar to Figure [1] 
we further mirror the flux maxima of SEPs on the right- 
hand half to the left-hand half with a blue dotted curve. We 
can obviously observe that with the same longitude sepa¬ 
rations between the solar source centers and the magnetic 
footpoint of the observer, the flux maxima of SEP events 
originating from solar sources located at eastern side of the 
observer footpoint are systematically larger than those of 



Figure 2: The flux maxima of the SEP events originating 
from 18 solar sources with different longitudinal locations 
relative to the magnetic footpoint of the observer at 1 AU, 
90° colatitude. With the same longitude separations be¬ 
tween the solar source centers and the magnetic footpoint 
of the observer, the flux maxima of the SEP events origi¬ 
nating from solar sources located at eastern (left) side of 
the observer footpoint are systematically larger than those 
of the SEP events originating from sources located at west¬ 
ern (right) side. 


the SEP events originating from sources located at western 
side. 


5 Relationship between Simulations and 
Observations 

Generally, a SEP event with a higher peak flux will reveal 
a larger intensity relative to other events through the en¬ 
tire evolution process including the rising phase and the 
decay phase (see the results, such as Figure 2 and Figure 
3, in 0). Generally, the probability of the SEP events be¬ 
ing observed by spacecraft near Earth’s orbit depends on 
the intensities of the particles in the events. An intense 
SEP event with high flux will experience relatively weak 
influences by the interplanetary structures (such as mag¬ 
netic cloud, heliospheric current sheet, corotating interac¬ 
tion region, etc.) during the propagation process in the he¬ 
liosphere. However, the not so intense SEP events with 
small fluxes will be significantly affected by the interplan¬ 
etary structures and then to weaken or even dissipate be¬ 
fore entering Earth’s orbit. Consequently, the intense SEP 
events would have larger probability of reaching the Earth 
and being observed by spacecraft at 1 AU; whereas the rel¬ 
atively weak SEP events would have lower probability to 
arrive at the Earth and to be recorded by spacecraft. 

As shown in Figure |2] the SEP events originating from 
solar sources located at eastern side of the magnetic foot¬ 
point of observer reveal larger flux maxima (corresponding 
to entire fluxes) than those originating from solar sources 
located at western side. Therefore, the SEP events origi¬ 
nating from solar sources located at eastern relative lon¬ 
gitudes would have higher probability to reach the Earth 
and to be observed by spacecraft near 1 AU. On the con¬ 
trary, the SEP events from solar sources on the west side 
relative to the magnetic footpoint of observer will be more 
easily influenced by the interplanetary structures and then 
to weaken or even dissipate during their propagation pro¬ 
cesses. As a result, the spacecraft at 1 AU would have a 
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larger probability to miss the SEP events originating from 
solar sources located at western relative longitudes. Ac¬ 
cordingly, taking into account a large data set with 78 so¬ 
lar proton events accumulated over a long time period, we 
can expect that with the same relatively longitudinal dis¬ 
tance, the number of the solar proton events originating 
from solar sources located at eastern side of the magnetic 
footpoint of observer would be larger than that of the solar 
proton events originating from sources located at western 
side, just as shown in Figure [T] Therefore, our numerical 
simulations explain and confirm the observational results; 
and the combination of simulations and spacecraft observa¬ 
tions demonstrates that there exists a longitudinally asym¬ 
metric distribution of solar sources of the proton events ob¬ 
served near 1 AU. 

6 Summary and Discussion 

In this work, we statistically investigate a large data set 
with 78 solar proton events accumulated over a long and 
uninterrupted time period 1996-2011. We find that with 
the same longitude separation between magnetic footpoint 
of observer and solar sources, the number of the proton 
events originating from the solar sources located at east¬ 
ern side of the magnetic footpoint of observer is systemat¬ 
ically larger than that of the proton events associated with 
the sources located at western side. Our simulation results 
show that with the same longitude separation between the 
solar sources and the magnetic footpoint of the observer, 
the flux of SEPs released from the solar source located east 
is systematically higher than that of SEPs originating from 
the source located west. Consequently, the SEP events orig¬ 
inating from solar sources located at the eastern relative 
longitudes would be detected and recorded with somewhat 
higher probabilities by the observer near 1 AU. Therefore, 
our numerical simulations testify and confirm the statisti¬ 
cal results of the solar proton events accumulated over 16 
years. 

Generally, as viewed from above the north pole of the 
Sun, the Parker spiral magnetic field lines originating from 
the solar surface would meander clockwise to somewhere 
with more eastern longitudes in the interplanetary space. 
The scenario can be seen in Figure^ As we can see, there 
is an east-west azimuthal asymmetry in the topology of 
the Parker interplanetary magnetic field. As a result, with 
the same longitude separation between the magnetic foot¬ 
point of the observer and the solar sources, the virtual dis¬ 
tance OA between the observer (O) in the interplanetary 
space and the solar source (A) on the left is shorter than 
that OB between the observer and the solar source (B) on 
the right. According to the theoretical and observational in¬ 
vestigations (e.g., Q), for physical conditions representa¬ 
tive of the solar wind, the perpendicular diffusion coeffi¬ 
cient would be a few percent of the parallel diffusion coef¬ 
ficient in magnitude. In general, SEPs would mainly trans¬ 
port along the average interplanetary magnetic field after 
they have been produced from the surface of the Sun or a 
CME-driven shock. Therefore, when the observer is con¬ 
nected directly to the solar source by the interplanetary 
magnetic field lines, the omnidirectional flux and peak flux 
of SEPs are larger than those in the cases where the ob¬ 
server is not connected directly to the solar sources. The 
farther the solar source center is away from the magnetic 
footpoint of the observer, the smaller are the omnidirec¬ 
tional flux and peak flux of SEPs observed. In addition 



Figure 3: Sketch to show the east-west azimuthal asymme¬ 
try in the topology of the Parker interplanetary magnetic 
field. With the same longitudinal separation between the 
magnetic footpoint of the observer and the solar sources, 
the virtual distance OA between the observer (O) in the in¬ 
terplanetary space and the solar source (A) on the left is 
shorter than that OB between the observer and the solar 
source (B) on the right. On account of OA < OB, the SEPs 
originating from the solar sources on the east (left) side rel¬ 
ative to the magnetic footpoint of the observer would have 
a lower probability to weaken or dissipate before arriving 
at the observer. 

to the parallel diffusion, the SEPs will experience perpen¬ 
dicular diffusion to cross the heliospheric magnetic field 
lines during their transport processes in the interplanetary 
space. This is why the SEP events can be observed by a 
spacecraft which is not directly connected to the accelera¬ 
tion regions. On account of the virtual distance relationship 
OA < OB, the SEPs originating from the solar sources lo¬ 
cated at eastern side of the magnetic footpoint of observer 
would have a lower probability to weaken or dissipate be¬ 
fore arriving at the observer. Consequently, with the same 
longitudinal distances between the solar sources and the 
observer footpoint, the fluxes of the SEP events associated 
with solar sources located east are systematically larger 
than those of the SEP events associated with sources lo¬ 
cated west. Therefore, the longitudinally asymmetric dis¬ 
tribution of solar proton events results from the east-west 
azimuthal asymmetry in the topology of the Parker inter¬ 
planetary magnetic field as well as the effects of perpendic¬ 
ular diffusion on the transport processes of SEPs in the he¬ 
liosphere. 
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